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ABSTRACT 

Aims. We report here on the first detection at X-ray wavelengths of the Supernova Remnant (SNR) G337. 8-0.1, carried out with the 

XMM-Newton Observatory. 

Methods. Using the X-ray observations, we studied the X-ray morphology of the remnant at diff'erent energy ranges, analysed the 

spectral properties and investigated a possible variable behavior. 

Results. The SNR shows a diff'use filled-center structure in the X-ray region with an absence of a compact source in its center. We find 

a high column density of A^h > 6.9x10^^ cm~^, which supports a relatively distant location (d >1 kpc). The X-ray spectrum exhibits 

emission lines, indicating that the X-ray emission has a thin thermal plasma origin, and is well represented by a non-equilibrium 

ionization (NEI) plasma model. The X-ray characteristics and well-known radio parameters show that G337. 8-0.1 belongs to the 

emerging class of mixed-morphology (MM) SNRs. 
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1, Introduction In this Letter we report the first detection at X-ray wave- 

lengths of SNR G337.8-0.1. The object is characterized by an 
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^ G337. 8-0.1 (also known as Kes 41) was discovered with the apparent centrally brightened X-ray morphology and a spectrum 
^ Molonglo telescope at 408 MHz (Shaver & Goss 1970). This ob- which suggests that the bright central emission is thermal in na- 
;^ ject, located in the southern galactic plane, was reported in the ture. The structure of the paper is as follows: in Sec. [2] we de- 
^^ Molonglo Observatory Synthesis Telescope (MOST) Catalog scribe the XMM-A^^w^o/z observations and data reduction. X-ray 
K> ■ (Whiteoak & Green 1996) as a bright small-diameter SNR and analysis and results are presented in Sec. [3 In Sec. H we discuss 
^ . distorted shell. The remnant is 6' x9' in size, distinctly elongated the implications of our results and summarize the main condu- 
ct in the northeast- southwest direction. Based on HI observations, sions. 
Caswell et al. ( 1 975) placed the SNR beyond the tangent point at 
7.9 kpc. Koralesky et al. (1998) detected maser emission in the 

f'^^tu^ -/5 ^"^ r'' ™P.iyi"g ^ far kinematic distance ofi 2.4 3. Observations and data reduction 

kpc. This fact is characteristic of a shock interaction with dense 

molecular gas (Green et al. 1997). The MOST flux density of The SNR G337.8-0.1 was observed on September 2005 by the 

S843MHz=18 Jy, combined with S408MHz=26 Jy (Shaver & Goss XMM-Newton X-ray satellite (Obs-Id. 0303100101). The obser- 

1970), gives a spectral index of a ~ -0.51 for the remnant. vation was centered towards the source (aj2ooo.o= 16^39"^06'.5, 

At X-ray frequencies the object was never detected by early ^J2000.o = - 46°57^58rO), and images were acquired withjhe 

X-ray missions. This non-detection is likely a result of the ab- EPIC MOS (iTurneij |2001|) and EPIC PN ( Struder et al, 2^3 

sorption of the typically soft (i.e., < 2 keV) X-ray emission by cameras. The observation was made through a ''thin'' filter, and 

the high column density of gas and dust in the Galactic plane. As in the full frame (FF) imaging mode. The observation was re- 

a part of an efl'ort to identify and study the X-ray emission of the trieved from the XMM-Newton Science Archive (XS AU and 

SNR G337.8-0.1 we have used the greatly enhanced X-ray sen- raw EPIC data were calibrated using the last version of the 

sitivity provided by the XMM-Newton telescope. An instrument Standard Analysis System (SAS-V7. 1 .2B To construct images, 

with these characteristics ofl'ers a unique opportunity to study spectra and light curves, we selected events with flag==0, and 

heavily obscured SNRs. patterns< 12 and 4 for the MOS and PN cameras, respectively. 



http://xmm.vilspa.esa.es/xsa/ 
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Fig. 1. XMM-Newton image of Kes 41 in the 0.5-10 keV energy 
range. The radio map observed with the MOST at 843 MHz is 
overlapped in yellow contours (Whiteoak & Green, 1996). Radio 
contours are indicated in steps of 0.2, 0.29, 0.38, 0.47, 0.56, 0.65, 
0.74 and 0.83 mJy beam"^ The resolution is ~ 43"x43" cosec 
(S). The red ellipse shows the region from which the X-ray spec- 
tra was extracted (see Sec. 3. 2). 



We derive Good Time Intervals (GTI) by the accumulation 
of background light-curves in the 10-15 keV energy band, which 
leads to a reduction of -15% in the net exposure time of the ob- 
servation. The net exposure time of the observation for all EPIC 
cameras finally reaches about 44.8 ks. Unfortunately, the object 
was negatively aff'ected by bad columns and CCD gaps in the 
EPIC-PN camera, which makes the spatial and spectral study of 
the X-ray emission from the SNR unreliable. In order to avoid 
systematics that aff'ect the reliability of our analysis, we made 
use of only the EPIC-MOSl/2 data. 



3. X-ray analysis of G337.8-0.1 

3.1. Image 

We use the clean event files to generate MOSl and M0S2 im- 
ages in the energy band [0.5-10] keV with a spatial binning of 
6.53 arcsec per pixel. In order to increase the signal-to-noise 
(S/N) ratio, we use the emosaic SAS task to merge the two 
MOSl/2 images. The corresponding set of exposure maps for 
each camera was prepared to account for spatial quantum eflft- 
ciency and mirror vignetting by running the SAS task eexmap. 
Exposure vignetting corrections were performed by dividing the 
superposed count image by the corresponding superposed expo- 
sure maps. We adaptively smoothed this image to a S/N ratio of 
10 using the SAS task asmooth. 

Figure 1 shows the X-ray image of the SNR G337. 8-0.1 in 
the 0.5-10.0 keV energy band with the radio contours at 843 
MHz superposed. The image reveals diff'use X-ray emission with 
an apparent filled-center structure and the absence of a compact 
source in its center. The X-ray peak of the SNR G337. 8-0.1 is 
located at (aj2ooo.o = 16^38^55^.7, ^j2ooo.o = -46°58'32r4), and 
as can be seen, the X-ray emission region is smaller than the 
radio structure. In fact, the X-ray emission fills the interior of 
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Fig. 2. Composite XMM-Newton images of G337. 8-0.1 at low, 
median and high energy (0.5-2 keV, 2-4 keV, 4-6 keV, 6-10 keV) 
for the combined MOSl/2 cameras. Note the lack of soft X-ray 
photons in the 0.5-2 keV energy range, while mostly of the X- 
ray emission is clearly observed in the 2-4 keV image and the 
semi-circular structure in the 4-6 keV band. 



the radio shell and lies mainly in the southwest part of the radio 
shell. 

To examine the morphology of the remnant in more detail, 
we have generated narrow-band images in the energy ranges 
0.5-2 keV, 2-4 keV, 4-6 keV and 6-10 keV. In Figure 2 we 
show the images in the four narrow energy bands. The source 
is almost undetected below 2 keV. The distribution of the 
emission in the 2-6 keV range is quite similar to the broad- 
band image. The centroid of the hard X-ray emission lies in 
the SW region of the remnant. The absorbed X-ray fluxes 
for the bands are: Fo.5-2keV= 1-3x10"^^ ergs"^ cm"^, F2-4kev = 
1.6x10"^^ ergs-^cm-^ F4-6keV=7.35xl0-^4 ergs'^cm'^ and 
F6-iokeV =2.9x10"^"^ ergs"^cm"^, which correspond to 5%, 
58%, 26% and 11% of the total X-ray emission, respectively. 
As can be seen, mostly of the X-ray emission originates in the 
2-6 keV band. 



3.2. Spectral analysis 

For the spectral analysis we used MPS 1/ 2 data. It was performed 
using the xspec package (lArnaudi Il996l) . Since the statistics of 
the source are not complete enough to perform a spatially re- 
solved spectral analysis, we extracted X-ray photon events for 
the whole source by using an elliptical region with a major axis 
of 5 arcmin and a minor axis of 3.2 arcmin, as indicated in 
Fig. 1 . The background region was taken from a nearby blank re- 
gion in the neighborhood of the source. Ancillary response files 
(ARFs) and redistribution matrix files (RMFs) were calculated. 
The spectra were grouped with a minimum of 12 counts per bin. 
The background- subtracted spectra of the MOS data are shown 
in Figure 3. At high X-ray energies (above 6 keV), the spectra 
show features with low statistical significance (~ 1 to 1.5 sigma), 
and are probably related to fluoresc ence lines in the backgroun d 
spectrum of the XMM-A^^w^6>^ (e.g. lDe Luca & Molendi[l2004 . 
Our analysis of the XMM-Newton EPIC spectra was per- 
formed using a NEI model (Mazzotta et al. 1998) (see Fig. 3). 
The model was afl'ected by a multiplicative absorption wab s 
(wABSisM + wabssnr) modcl (iMorrison & McCammonl Il983h . 
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channel energy (keV) 

Fig. 3. The XMM-Newton MOSl/2 spectrum of the diffuse X- 
ray emission of Kes 41. The solid line indicates the best fit for 
the NEI model (see Table 1). The soft X-ray excess are proba- 
bly due to foreground emission. The position of the most intense 
transition lines of Sxv and Sixiv elements are indicated. The to- 
tal number of photons is 3725 and 2539 for MOSl and M0S2, 
respectively. 

Table 1. X-ray spectral parameters of the SNR G337.8-0.1 



Model: 


NEI 


Parameter 


Value 


Icr error 


N!?^^ [cm-'] 
N^^ [cm-2] 


2.2x10^' 


fixed 


4.7x10^2 


1.0x10^2 


kT [keV] 


1.4 


0.2 


Abundance^ 


0.6 


0.3 


log(T) [scm-3] 


11.0 


0.3 


Normalization 


3.4x10-3 


1.1x10-3 


Flux [ergcm"^ s"^] 


7.5x10-12 


1.8x10-12 


X^/y-d.o.f. 


1.08/380 


— 



Notes: The ISM Hi absorption column density (N™) was com- 
puted by following Dickey & Lockman (1990) HI in the Galaxy, and 
kept fixed in the fit procedure. The logarithm of ionization timescale 
log(T) is also written as log(net). Flux is absorption-corrected in the 
0.5-10.0 keV energy range. Normalization was calculated accord- 
ing to: 10~^^/47ir^ J nenndV. t The abundance was adopted from 
I Anders & Grevessd (.1989.) and left as a free parameter in all our fits. 



The goodness of the model fit was derived according to the;^^- 
test statistics. C-statistics was also applied to the spectral fit, giv- 
ing similar results. Results of the spectral analysis are summa- 
rized in Table 1. The presence of excess soft emission is not 
surprising, given that Kes 41 is located in a crowded region of 
the galaxy at a distance of between 7 to 12 kpc. Emission-line 
features corresponding to the blend of S xv and Si xiv are clearly 
visible at ~ 2.4 keV. Possible detection of others elements like 
Arxvii and S xvi emission lines appear blended around 3.1 keV. 
Although the presence of these lines support a thermal nature for 
the X-ray emission, the absence of significant Fexxv at 6.7 keV 
does not allow us to discard a weak non-thermal contribution at 
high X-ray energies. 

To obtain a statistical assessment of the X-ray variability of 
the SNR G337.8-0.1, we use the 44.8 ks EPIC-PN observation 
to compare the time arrival distribution of source photons by 
means of the Kolmogorov-Smirnov (KS) test (Press et al., 1992). 
We use an extraction region centered on the whole X-ray emis- 
sion. We see no significant pulsed signal with a period greater 



than twice the read-out time of the EPIC-PN camera in the FF 
mode (73.3 ms), which corresponds to a Nyquist limit of 0.146 
s. At first glance, the absence of any significant X-ray variability 
and a thermal origin for the X-ray emission allows us to discard 
an energetic pulsar associated with the SNR. 

4. Discussion 

The SNR G3 37. 8 -0.1 is characterized by an apparent centrally 
filled X-ray morphology not well correlated with the radio shell, 
and thermally dominated X-ray emission. Thus, Kes 41 could 
be classified as a thermal composite (or mixed-morphology) 
remnant. Moreover, at radio frequencies, the emission is non- 
thermal (Whiteoak & Green 1996) and can be interpreted as 
synchrotron radiation from accelerated high-energy electrons at 
the SNR shock. These characteristics place the remnant in the 
emerging class of SNR that includes objects such as W44 (Jones 
et al. 1993), 3C 391 (Rho & Petre 1996); G272.2-3.2 (Harrus et 
al. 2001) and G290.1-0.8 (Slane et al. 2002). A fist containing 
SNRs with these characteristics was proposed by Rho & Petre 
(1998). 

Although the evolutionary states that led to the observed 
properties in these SNRs is not well understood, it is widely 
believed that these peculiar characteristics are linked to inho- 
mogeneities in the ISM. Several possible scenarios have been 
introduced in the past to explain thermal X-ray radiation inside 
radio shells of SNRs, i) cloudlet evaporation in the SNR inte- 
rior (White & Long 1991), ii) thermal conduction smoothing out 
the temperature gradient across the SNR and enhancing the cen- 
tral density (Cox et al. 1999), iii) a radiatively cooled rim with a 
hot interior (Harrus et al. 1997), and iv) possible collisions with 
molecular clouds (Safi-Harb et al. 2005). 

G337. 8-0.1 lies in a complex region of the sky, possibly 
adjacent to a very massive giant molecular cloud (Dame et al. 
1986) and close to the associated maser emission detected at -45 
km s"^ by Koralesky et al. (1998). If the SNR is indeed propa- 
gating in a region with a large number of clouds, the most suit- 
able scenario for describing the SNR evolution is that developed 
by White & Long (1991). In such remnants, the clouds evapo- 
rate by saturated conduction (Cowie & McKee, 1977) and the 
X-ray morphologies may be very diff'erent from ordinary shell- 
like SNRs. 

Based on the X-ray properties of G337.8-0.1, specifically 
the filled-center morphology and detected lines in the spectrum, 
a thermal interpretation of the X-ray emission is the most plau- 
sible to describe the SNR evolution. Using the X-ray image of 
G337. 8-0.1 we estimated the volume V of the X-ray emitting 
plasma, assuming that the plasma fills an ellipsoid similar to the 
region from which the X-ray spectra was extracted (with a diam- 
eter of 5'x3'.2). For a distance range from 7 to 12 kpc, we then 
obtained a volume of (2.8-14.8)xl0^^ cm^. Based on the emis- 
sion measure (EM) determined by the spectral fitting, we can 
estimate the electron density of the plasma, n^, by ne= ^JEM/V, 
which varies between 0.20 cm"^ and 0.26 cm"^. In this case, 
the number density of the nucleons was simply assumed to be 
the same as that of electrons. The range of age t was then de- 
termined from the ionization timescale, r, by t=T/ne. Therefore, 
the elapsed time after the plasma was heated is within (12.000- 
16.000) yr. The total mass of the plasma Mtotai was estimated by 
Mtotai=^e^ niH = (5 - 32) Mq, where mn is the mass of a hydro- 
gen atom. These results are consistent with the expected values 
for a middle-age SNR. 

One of the most detailed studies of middle-age and center- 
filled thermal emission in a SNR interacting with a dense molec- 



Combi et al: XMM-Newton observations of G337. 8-0.1 



ular cloud is the W44 model (Cox et al. 1999; Shelton et al, 
1999; Shelton et al. 2004). The proximity to molecular clouds 
of this SNR suggests that the environmental interstellar mat- 
ter, where the SNR is propagating, is relatively dense. Using 
Cox's model we have roughly estimated some dynamical char- 
acteristics (i.e., the shock radius R^ and shell velocity Vsheii) for 
G337.8-0.1. Assuming an energy for the explosion of ~ 10^^ 
ergs, we found that 7?s= (8-10) pc and Vsheii= (180-200) km s"^. 
These values are consistent with the angular size and range of 
distance for the SNR under consideration. Thus, G337. 8-0.1 
could be, a SNR with characteristics similar to W44. 

Finally, we have compared the X-ray luminosity of 
G337. 8-0.1 with other proto-typical mixed-morphology SNRs. 
The unabsorbed X-ray luminosities in the 0.5-10 keV band of 
G337. 8-0.1 for a range of distance between 7 and 12 kpc is 
Lx ~ 4.5x10^"^ erg/s and 1.3x10^^ erg/s, respectively. The first 
value is similar to the X-ray luminosity of W28 (Lx~6xl0^'^ erg 
s-^ Rho & Borkowski, 2002), G272.2-3.2 (Lx ~ 5.2x10^^ ^^^ 
s"^ ; Harrus et al. 2001), and Kes 32 (Lx ~ 3.2x10^^ ^^^ ^-i . yj^k 
2008) for the same energy range. This suggests that a distance 
around 7 kpc for G337. 8-0.1 is most probable. 

In summary, in this work we reported the first detection at 
X-ray wavelengths of the SNR G337.8-0.1. The radio and X- 
ray properties show that the object displays the three basics at- 
tributes found in mixed-morphology SNRs, i.e. centrally filled 
X-ray morphology with lines in its spectra, a shell-like appear- 
ance at radio frequencies, and the absence of a prominent, cen- 
tral, compact source in radio and X-ray energies. If the SNR is 
indeed propagating in a cloudy medium, the interaction of the 
shock- front with the adjacent molecular material could be re- 
sponsible for the gamma-ray emission detected by the EGRET 
telescope. 

The detailed properties of the SNR are still unknown because 
of the low photon statistics in the XMM-Newton data. Future X- 
ray observations with a higher spatial resolution and good pho- 
ton statistics like Chandra or Suzaku will reveal the relation be- 
tween the X-ray intensity peak and the radio emission. 
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